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Abstract 



^^ I Dimensional regularization is applied to the Lippmann-Schwinger equation 

for a separable potential which gives rise to logarithmic singularities in the 
^ ' Born series. For this potential a subtraction at a fixed energy can be used to 

renormalize the amplitude and produce a finite solution to the integral equa- 

r> I tion for all energies. This can be done either algebraically or numerically. 



In the latter case dimensional regularization can be implemented by solving 
the integral equation in a lower number of dimensions, fixing the potential 
strength, and computing the phase shifts, while taking the limit as the num- 
ber of dimensions approaches three. We demonstrate that these steps can be 
carried out in a numerically stable way, and show that the results thereby 
obtained agree with those found when the renormalization is performed alge- 
braically to four significant figures. 



I. INTRODUCTION 

One difficulty for standard treatments of hadronic reactions is that form factors are intro- 
duced at hadronic vertices in order to regulate integrals which would otherwise be divergent. 
This procedure reflects the substructure of hadrons which gives them a finite extent, and 
hence a form factor. However, any field theory upon which such calculations are based 
will necessarily be non-local. The implementation of basic field theoretic principles, such as 
causality and electromagnetic gauge invariance, is quite involved in such field theories [|T],3. 
Methods which impose gauge invariance on an amplitude containing hadronic form factors 
have been formulated, ^-^ but they are intrinsically non-unique, since they constrain only 
the longitudinal part of the photon's coupling to the hadronic system. 

Some of these difficulties can be resolved by using dimensional regularization (DR) |§-|^ 
to render divergent integrals finite. This is the method of choice for dealing with the infinities 
which arise in perturbative field theoretic calculations. However, there are very few studies of 
the application of dimensional regularization to integral equations. One notable exception is 
the recent application of DR to a Schwinger-Dyson equation in quenched QED by Schreiber 
et al. ||10|. In this paper we adapt their method and apply it to the Lippmann-Schwinger 
(LS) equation. 

We consider one member of the class of potentials in which one subtraction at a fixed 
energy results in a finite integral equation at all energies. The potential we choose is separable 
and so provides a toy problem where the renormalized amplitude can be derived algebraically 
in a straightforward way. The amplitude thereby obtained can then be compared with that 
found by the numerical solution of the dimensionally-regulated integral equation. We will 
show that with careful numerical work the phase shifts obtained via these two different 
methods are in excellent agreement with each other. This proves that many of the difficulties 
associated with implementing dimensional regularization numerically can be overcome, and 
hence that such "numerical dimensional regularization and renormalization" facilitates the 
extraction of finite phase shifts from a particular class of divergent potentials. 



In Sec. H we outline the general framework for implementing dimensional regulariza- 
tion and renormalization techniques numerically in the LS equation. We show how the LS 
equation can be analytically continued to lower numbers of dimensions, D = 3 — e, thus 
ameliorating the divergences which appear in three dimensions. The application of a renor- 
malization condition at one energy in 3 — e dimensions can then lead to an amplitude which is 
finite at all energies in this lower number of dimensions, and thence to an amplitude which is 



well-defined as the limit e ^ is taken. In Sec. |ITl| the special case of a separable potential, 
which leads to a divergent amplitude as e — i> is considered. We obtain the amplitude for 
this potential for e > 0. Renormalizing the amplitude by demanding that the binding energy 
of the deuteron be reproduced leads to a divergence in the inverse of the "bare" coupling 
which appears in the original separable potential. We show that in DR this divergence goes 
as 1/e. We then observe that for this particular potential, it is not necessary to use DR, 
since the problem can be subtractively renormalized. That is to say, a subtraction at one 
energy, which can be carried out algebraically, renders the amplitude finite at all energies. 

However, it is not generally true that such algebraic renormalization of an integral equa- 
tion is possible. Therefore, in Section ^V\ the general renormalization procedure of Section |I| 
is applied directly to the LS equation in 3 — e dimensions. In this approach no subtraction 
is carried out, and no assumption is made about the form of the potential. First, the LS 
equation is solved numerically in D = 3 — e dimensions. Then, demanding a specific value 
of the amplitude at one energy implies a value of the bare coupling in D dimensions. The 
amplitude can then be computed numerically for this value of the coupling at arbitrary 
energy and physical quantities, such as phase shifts, extracted in D dimensions. This must 
be done over a range of values e and then the limit e ^ can be taken. Finally, in Sec. 
we present some concluding remarks. 



II. THEORETICAL ANALYSIS 

Dimensional regularization is often used to identify the infinities encountered in pertur- 
bative calculations 0-|^. The procedure involves replacing integrals that are infinite in n 
dimensions by the same integrals in a lower dimension n — e, where they have a finite result. 
The function defining this integral, say /(e), can then be analytically continued back to n 
dimensions. In this way the logarithmic singularities of / can be isolated and then renormal- 
ized by adjusting the bare parameters of the theory. This is the method of choice for dealing 
with the infinities which arise in perturbative quantum field theory calculations, but it has 
not been widely applied to non-perturbative integral equations that contain divergences. 



It was used in the Lippmann-Schwinger equation for a divergent potential in Refs. ||ll| , p!2 
but in the cases considered there the amplitude can be derived analytically. In the case 
of a general divergent potential, this will not be possible, and DR must be implemented 



numerically. Recently, Schreiber et al. |T^ have demonstrated the feasibility of such a pro- 
cedure, successfully employing DR in the numerical solution of a divergent four- dimensional 
Schwinger-Dyson equation for the electron mass in quenched QED. In this section we explain 
how to adapt the methods of Schreiber et al. for use in the three-dimensional Lippmann- 
Schwinger equation. 

Given a potential V the momentum-space LS equation for the scattering amplitude T 
for two particles of mass m is: 

r(k, k'; E) = r (k, k') + / d'k" V(k, k") J r(k", k; E) , (1) 

J E — k'^^/m + IT] 

where 77 is a positive infinitesimal. To dimensionally-regulate this equation we must continue 
all quantities in it into 3 — e dimensions. If the potential V is central, and depends only on 
the angle between k and k', then the operator L^ commutes with the Hamiltonian. It follows 
that V and T can be expanded using the eigenfunctions of the rotation operator in 3 — e 
dimensions, thus reducing Eq. (|I]) to a one- dimensional integral equation. This expansion 
can be quite involved (see, for instance, the analogous expansion in 4 — e dimensions of 
Ref. [0). 



Here we follow a slightly different approach. Although the Lippmann-Schwinger equation 
is nominally a three-dimensional integral equation, in most instances it is written and solved 
as a one-dimensional integral equation for a given partial wave L In other words, making 
the expansion: 

T(k, k'; E) = ^ 5:(2£ + l)T\k, k'- E)Pi{k ■ k') , (2) 

and an equivalent expansion for V , leads, as is well-known, to the one-dimensional integral 
equation: 

T\k, k'- E) = v\k, k') + J dk" k"' ^}!:,f\. nk", k'- E) . (3) 

To dimensionally-regulate this equation is much easier, since the angular part of the 
problem has already been dealt with in the original number of dimensions, in this case, 
three. It only involves replacing the measure appropriate to three dimensions, by one for 
(3 — e) dimensions, i.e. 

Tf (A:, k'- E) = V\k, k') + 7 dk" k"^'-^^ ^ K^^fl ■ T'Ak". k'; E) . (4) 

J L — k'^^/m + IT] 

It is worth pointing out here that if the partial-wave expansion (^ is carried out in 3 — e 
dimensions, rather than in three dimensions, then the potential V^ ends up depending on e 
as well as i. (See, for instance, Ref. 0].) However, if we are only concerned with isolating 
the divergences of the integral equation as e — ;> 0, then we are free to take the limit as e ^ 
first in V^, since that limit is non-singular. Thus, choosing to dimensionally regularize the 
partial-wave-expanded equation, rather than the full three-dimensional equation, may be 
unfamiliar to those used to DR in perturbative settings, but it is justified. 

Now consider a potential V^{k', k) which leads to divergences in the Born series for 
Eq. (^). The integral equation (^ will yield finite answers for the terms in that series, 
provided that e is taken large enough. Equation (^) has been regularized. The next step is to 
renormalize it. To that end we pick one of the parameters of the potential, e.g. its strength, 
which we denote by A, and regard it as a function of e, A = A^. A renormalization condition 



is then chosen. This can be represented as the demand that the amphtude, T^{k, k'; E), have 
some specified value at a particular kinematic point. For instance, we might demand that 
the amplitude has a pole aX E = ~B, so that a bound state exists at that energy. Such a 
renormalization condition implicitly defines the function A^. The amplitude T^{k', k; E) can 
then be computed at the on-shell point: k'"^ = k"^ = mE, for a variety of energies, with this 
particular choice of A^. This will certainly give finite results for e > 0, provided that e is big 
enough. Furthermore, if it is enough to renormalize only the strength of the potential, such 
a calculation will produce finite values for the on-shell amplitude and hence for the phase 
shift even in the limit e — > 0. This "numerical renormalization" program has been carried 



out before in the Lippmann-Schwinger equation with other regulators (see, e.g. |jT3|-0), but 
here we will implement it for DR. 

In the next section we illustrate how this numerical renormalization proceeds in the case 
of a simple separable potential. We show that the strength of the potential can always be 
adjusted so that the potential supports a bound state with energy E = —B. This strength 
cancels the divergence in the Born series integrals, and so a finite scattering amplitude can 
be extracted. 

III. THE SEPARABLE POTENTIAL AND ALGEBRAIC RENORMALIZATION 

Consider a separable potential with a form factor that does not guarantee the convergence 
of the integrals needed to determine the scattering amplitude. An example of such a potential 
is (from this point on we drop the angular momentum label t) 

V{k,k')=g{k)\g{k') with g{k) = j^^^^^^ . (5) 

The corresponding off-shell T-matrix is given by 

T{k,k'-E^)=g{k)T{E+)g{k'), (6) 

with 



HE)]-' = X-' - I{E) , (7) 

where the integral I{E) is 

I(E) = m[dke—^^ . (8) 

^ ' J mE-k^ + in ^ ^ 

' 

Since g{k) — > k^^^"^ as k ^ oo, the integrand in Eq. (^ goes as k^^ for k -^ oo, and the 

integral I{E) has a logarithmic divergence. 

Of course, if we had inserted the potential (0) in the dimensionally-regulated Lippmann- 

Schwinger equation, then the only thing that would have changed in the above analysis is 

that equations (0) and (H) would have become: 

H{E)]-' = X;'- UE) (9) 

h{E) =mfdk k'-^ y^''^}' . . (10) 

^ ^ J mE -k^ + 171 ^ ^ 

' 

This integral is now finite for any e greater than zero. The logarithmically divergent part 
of Jg will stand revealed as a 1/e pole. Since A^ is a function of e here this pole can be cancelled 
by making an appropriate choice of A^. Specifically, here we achieve this renormalization by 
demanding that one observable be reproduced, i.e. that the ^Si A^A^ channel have a bound 
state at energy Eb = —2.2246 MeV. This is done so that our interaction can be regarded as a 
toy model of the A^A^ potential in the deuteron channel. This model is extremely simple, but 
it has the only two features we are concerned with here: it leads to logarithmic divergences 
in the Born series, and it has a shallow bound state. The renormalization condition which 
fixes the position of this bound state may be expressed as follows: 

[t,{E)]-' = \-' - h{E) ^ Q as E^Eb = -B <0 . (11) 

To see how this removes the divergent part of Ie{E), we write the form factor g{k) as 

g[k) = —j= h{k) with h{k) — > 1 as k ^ oo . (12) 

Our singular integral can then be written as 



h{E) = -m J dk /^(i-) ,, ^^^^y . (13) 

^ ^ J k^ -mE-iri ^ ^ 

' 

We now define an integral Ie{E) that retains the singular part of Ie{E) in the hmit as e — ;► 0: 

W - -m I * ,,_„^_^ . (14) 

' 

This allows us to write Ie{E) in terms of a subtraction as 

IHkW 1 



I,{E) = -m I dkk^^-'^ I -J 

n V 



+ hi-B) 



mE — irj k"^ + mB 
= J,{E) + h{-B), (15) 






where the integral Je{E) is well defined in the limit as e -^ 0. Making use of the identity: 

r^ _i r(^)r(a-^ 

(r2 + A)° 2 A"-(/3+i)/2 r(a) 
we have 

f. m m r (1 - 6/2) r (1 + 6/2) m 

Ie{-B) = , > as 6^0. (17) 

e [mBy'^ e 

This isolates the singularity of Ie{E) in the limit 6^0. We now need to renormalize the 
strength of the potential A^ in order to get a finite amplitude. 
To do this we write the condition for a bound state (0) as 

K' = U-B) (18) 



I.e. 

which makes it clear that A^^ has both a divergent part, proportional to 1/6 in DR, and 
other, finite, parts, that are dependent upon the particular bound-state energy we are trying 
to reproduce. Equation (|I9|) gives A^ for a given e and from that we can then calculate Te{E) 
and the T-matrix. 

Alternatively, we can use Eq. ([T8| ) for A^ to write: 
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Tr\E) = U-B) - h{E) . (20) 

Inserting Eq. (|10|) for I^ and combining the integrands we see that: 

which imphes that: 



7;(fc',A;;i?)=^(fc')|^^(fc), (22) 



where: 



is a finite quantity. Since S^{E) is finite, the hmit as e goes to zero can be taken smoothly, 
and we never even need to consider any quantity other than So{E), when the regulator has 
been removed altogether. It is now trivial to obtain phase shifts from Eqs. (|22| ) and (p3[). 
Their relationship to the on-shell amplitude is just: 

^ ' Real T{E; ko, ko) ^ ^ 

where k^ = mE is the on-shell momentum. 



IV. NUMERICAL RENORMALIZATION 

The question motivating this study is whether numerical techniques can be found that 
give the amplitude T^{k', k; E) numerically when the LS equation is solved in less than three 
dimensions. These techniques must be stable enough to allow the limit e ^ to be taken. 
For the simple potential (^) the result of any such procedure must be given by Eqs. ( p2D and 
(p3|), where an "algebraic renormalization" making explicit use of the potential's separability 
has been carried out. 

In this section we implement an alternative strategy to this algebraic renormalization, 
that of Section |T|, which we call "numerical renormalization" . We solve the partial- wave 
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LS equation in (3 — e) dimensions, Eq. (H), directly. The strength A^ is adjusted to give 
the correct bound-state energy, i.e. so that T^{k,k'] E) has a pole aX E = —B. Once this 
is done we calculate the phase shifts as e — ;► 0. Here the dimensionally-regularized integral 
equation is being used to renormalize the strength of the potential. We do not have to 
explicitly carry through the subtractive, algebraic renormalization discussed in Sec. |ITI|, and 
the dimensionally-regularized LS equation should give phase shifts as a function of e which 
are well-behaved in the limit e ^ 0. 

The two signals we look for to determine the success of our numerical implementation of 
DR for our divergent integral equation are: 

1. Af must have the behavior (|T9|) in the limit e — »• 0; 



2. The results for phase shifts from numerical renormalization must be stable in the limit 

These two tests should be applicable to any dimensionally-regularized integral equation. 
However, in the case under investigation here the results of Eqs. (^^ and (^) provide 
an additional check on the accuracy of this numerical renormalization using dimensional 
regularization. 

To perform such tests, we need to solve the integral equation for each value of e for both 
the bound state problem and the scattering problem with considerable accuracy. To achieve 
this accuracy, we first split the integral in the Lippmann-Schwinger equation arbitrarily into 
a piece from zero to some km and a piece from km to infinity. Then, on the interval [km-, oo) 
we make a change of variables to: 

*^.log(|-). (25) 

thereby enabling us to employ a logarithmic mesh which ensures the correct numerical 
integration when we obtain l/A^. The variable transformation (^) introduces a factor e~* 
into the integrand, and so Gauss-Laguerre quadratures are chosen for the integration from 
km to infinity, since they naturally build in this factor. 
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Meanwhile, ordinary Gauss-Legendre quadratures are employed for the integration on 
[0, km]- This was done in two ways when we performed "numerical renormalization" : 

1. We deformed the contour of integration from the real fc-axis to that depicted in Fig. |l|. 
Along the contour from zero to k^.e'^'^ we divide the interval into two parts with hq 
quadratures for —>■ 2kQe~^'^ and ni quadratures for 2kQe~^'^ -^ kme~^'^. Finally, we 
take n2 quadratures for the part of the contour that returns us to the real axis, and 
n^ for the interval [k^, cxd). In this way we can optimize the four different regions of 
the integration independently. For the determination of A by the requirement that 
the potential supports a bound state, we take no = 0, i.e. the contour of integration 
corresponds to taking k^ = 0. We have optimized the number of quadratures on each 
of the intervals of integration as well as the angle and the point at which the contour 
returns to the real k axis, km- We have found it necessary to take no = 16, rii = 
80, n2 = 10, na = 15, = 0.7 and km = 50 fm^^. We could have used a smaller 
number of quadratures on each interval, but to establish that DR is valid for the LS 
equation and for the two methods to give identical results, we have not economized on 
the number of quadratures. 

2. Instead of solving for the T-matrix we solved for the K- matrix and used the relation- 
ship: 

T\E+) = K\E) - - mm/to K\E) T\E+), (26) 

to relate the two on-shell, and extract the phase shifts. The K-matrix is, of course, 
a purely real quantity, but the integral equation defining it has a principal-value sin- 
gularity at k = kQ, the on-shell point. To deal with this we place uq Gauss-Legendre 
quadratures, distributed symmetrically about k = ko, on the interval [0, 2fco]- We then 
also place rii Gauss-Legendre quadratures on [2ko, km] and n^ Gauss-Laguerre quadra- 
tures on [km, C)o). Here we found it sufficient to take no = 16, rii = 60, na = 15, and 
km = 50 fm~^. 

11 



Each of these methods for solving the integral equation is accurate to four significant figures, 
provided that e > 10^^. For very small e only the second mesh produces results which are 
this stable. 

Note that only the first integration technique was employed to do the integration in 
Eqs. (^) and (PBD, but there the precise details are less important, since the integral in 
question is finite. 

The first check that the dimensional regularization of this integral equation is being done 
correctly is to see that the bare coupling A extracted by imposing the condition that there 
be a bound state at E = —B on Eq. (^ does behave according to Eq. (|^). In Fig. |] we 
plot the function l/A^ versus 1/e that is found when this condition is imposed on Eq. (^) in 
the case of the potential (1). We chose B = -2.225 MeV, and /? = 1.4489 fm^^ The slope 
of the curve, is indeed, m, the nucleon mass, as per Eq. (plQl). 

Second, we can examine the convergence of the phase shifts with e. In Figs. ^ and ^ we 
present the results obtained from the Lippmann-Schwinger equation (^ with the potential 
(|) for the phase shifts as a function of 1/e. In fact the points calculated using the two 
different meshes described above are indistinguishable from one another on the scale shown. 
The calculation has been done for the nucleon-nucleon system in the '^Si channel at labora- 
tory energies of 24 and 352 MeV. The strength of the potential Ae varies with e as displayed 
in Fig. H, having been adjusted so that the potential supports a bound state with a binding 
energy of 2.225 MeV. 

A detailed comparison of the phase shifts resulting from the solution of the dimensionally- 
regularized LS equation with those found using the algebraic results of the previous section 
shows that the agreement is good to four significant figures until we get to very small (< 10^^) 
values of e. Long before that the phase shifts have converged as a function of e. This 
comparison is presented in in Tables | and y, where the the phase shifts obtained by these 
two methods at three energies, as well as the bare coupling A^, are shown as a function of e. 
Table | shows the results found by solving the homogeneous, dimensionally-regularized LSE 
to get Ae and then using that A^ to calculate phase shifts in the dimensionally-regularized 
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LSE Q. Table |ll| gives the result for phase shifts from Eqs. (p2D and (^) as well as the 



result (|T8D for A^. We observe that there is agreement between the algebraic and numerical 
renormalization to four significant figures. Furthermore, there is convergence in the phase 
shift as e ^ to five significant figures, which is beyond the numerical accuracy of this 
calculation. All this indicates that the dimensionally-regularized integral equation is giving 
a unique solution. 

Finally, we should mention that there is some sensitivity to the value of km that is chosen. 
This is a numerical effect, and reflects the wide spacing of quadratures in the logarithmic 
mesh above k = km- In Figure |^ we plot 6 at E\s,h = 352 MeV for e = 6.1035 x 10^^ over a 
range of kmS, using meshes of the second type described above. Provided km is large enough 
the results are stable to four significant figures. The variations in other phase shifts, and in 
Ae, are smaller than that displayed in this plot. 

V. CONCLUSION 



From the above analysis and that of Ref. [^ we conclude that it is possible to use di- 
mensional regularization to render divergent integral equations with logarithmic divergences 
finite, provided that a renormalization condition is imposed in order to fix one of the param- 
eters of the potential. Although in the present investigation we have chosen the bound-state 
energy to fix the strength of the potential, we could just as easily have used the scattering 
length or the value of the on-shell amplitude at some finite energy as the renormalization 
condition. 

At this point one might ask whether the ideas discussed here can be profitably employed 
if power-law divergences are involved. Of course, such divergences do not appear explicitly 



when DR is implemented in analytic calculations (see, e.g., Ref. |1T1|). However, the numer- 



^The results displayed in Table | were obtained using the K-matrix method and the second of the 
two meshes described above. 
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ical techniques discussed above simply will not eliminate ultra-violet divergences of degree 
greater than zero. The reason for this lies in the way such divergences are eliminated in 
"standard" DR. There the offending integral is analytically continued into a region where 
e is large enough so that the divergences no longer appear. The resulting analytic form 
is then defined as the value of the integral in the region where the integral was formally 
divergent. It is this definition that eliminates the power-law divergences, and, in contrast 
to the case where logarithmic divergences are present, the limit e — ;> cannot be taken until 
this additional step is made. Since the work of this paper relies on being able to straight- 
forwardly take the limit e — >■ it is not clear that power-law divergences can be eliminated 
in "numerical" DR. It is possible that a numerical procedure analogous to the analytic one 
just described can be developed to make sense of integral equations in which power-law 
divergences appear in the scattering series. However, the working out of such a scheme is 
beyond the scope of this paper. 

Although the present analysis is restricted to a simple toy model which has only log- 
arithmic divergences, we see no reason why it cannot be easily extended to other integral 
equations of interest. Two examples which we believe to be particularly important are: 

1. Effective field theory treatments of neutron-deuteron scattering in the doublet channel. 
The leading-order effective field theory calculation produces Faddeev equations whose 
kernel does not go to zero fast enough as /c — >■ cxd to make the integral equation 
well-behaved. Normally this equation is regulated by a cutoff, but the techniques 
discussed here could also be used. This problem is particularly interesting because 
it has recently been shown that unless a three-body force is added to the leading- 
order effective field theory calculation the resultant amplitude is unduly sensitive to 
the value of the cutoff |T^. This three-body force introduces a new parameter into 
the calculation, which is fit to the NNN doublet scattering length, ^a. The integral 
equation is then properly renormalized, in the sense that its solutions are no longer 
sensitive to physics at short distances in the NNN system. However, formally it still 
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contains divergences. 



2. The Bethe-Salpeter equation for pion-nucleon scattering [|T8|,|T9|. In this case we would 
hope to use DR rather than introduce form factors when solving the integral equation. 
This might facilitate the introduction of electromagnetic couplings into the it — N 
scattering problem. 
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TABLES 



e 


A (fm) 


5(24) 


5(96) 


5(352) 


1.0 


0.11928 


57.352 


37.151 


22.553 


2.5 xlO"^ 


4.9817 xl0"2 


89.360 


67.566 


49.740 


6.25 xlO-2 


1.30340 xl0~2 


97.748 


76.120 


58.067 


1.5625 xlO-2 


3.27867 xlO-3 


99.855 


78.299 


60.224 


3.90625 xlO-3 


8.20662 xlO-4 


100.38 


78.848 


60.770 


9.76563 xlO"^ 


2.05224 xlO"^ 


100.51 


78.985 


60.906 


2.44141 xlO"^ 


5.13094 xlO"^ 


100.55 


79.019 


60.941 


6.10352 xlO-5 


1.28276 xlO-5 


100.56 


79.028 


60.949 


1.52588 xl0~5 


3.20691 xl0~6 


100.56 


79.030 


60.951 


3.8147 xlO-6 


8.01728 xlO-^ 


100.56 


79.031 


60.952 


9.53674 xlO'^ 


2.00432 xlO"^ 


100.56 


79.031 


60.952 



TABLE I. Couplings and phase shifts for various values of the distance e, away from three 
dimensions using the subtraction technique. All quantities are numerically accurate to four signif- 
icant figures. Generated with uq = 16; n-i = 60; n^ = 15, kmax = 50 fm^ . 



e 


A (fm) 


5(24) 


5(96) 


5(352) 


1.0 


0.119286 


57.356 


37.154 


22.555 


2.5 xlO"^ 


4.98162 xlO-2 


89.355 


67.563 


49.738 


6.25 xlO-2 


1.30340 xlO-2 


97.742 


76.114 


58.059 


1.5625 xl0"2 


3.27867 xlO-3 


99.849 


78.297 


60.222 


3.90625 xlO-3 


8.20662 xlO^*^ 


100.38 


78.845 


60.768 


9.76563 xlO-4 


2.05223 xlO'^ 


100.51 


78.983 


60.905 


2.44141 xlO"^ 


5.13094 xlO-5 


100.54 


79.017 


60.939 


6.10352 xlO-5 


1.28276 xlO-5 


100.55 


79.026 


60.948 


1.52588 xlO-5 


3.20691 xlO-6 


100.55 


79.028 


60.950 


3.8147 xl0^6 


8.01728 xlO"'^ 


100.55 


79.028 


60.950 


9.5367 xlO-^ 


2.00432 xlO-^ 


100.55 


79.028 


60.950 



TABLE II. Couplings and phase shifts for various values of the distance e, away from three 
dimensions using the algebraic solution of the separable potential and contour rotation. All quanti- 
ties are numerically accurate to five significant figures. Generated with no = 16, rii = 80, n2 = 10, 
713 = 15, 4> = 0.7 and km = 50 fm^"*^. 
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FIGURES 
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FIG. 1. The rotated contour of integration. Here k^ is the on-shell momentum 



(1/^J(fm') 




FIG. 2. The behaviour of the inverse of the bare couphng A as a function of 1/e, as obtained 
using the dimensionally-regularized Lippmann-Schwinger equation 
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FIG. 3. Phase shifts at EL^b = 24.0 MeV for the solution of the LS equation as a function of e 
for e ^ 0. 
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FIG. 4. Phase shifts at ELab = 352.0 MeV for the solution of the LS equation as a function of 
e for e ^ 0. 
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FIG. 5. Phase shifts at ELab = 352.0 MeV with e = 6.1035 x 10"^ for the solution of the LS 
equation as a function of km- 
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